Anisotropic Diffusion-Limited 
Reactions with Coagulation and Annihilation 



Vladimir Privman, Antonio M. R. Cadilhe, and M. Lawrence Glasser 

Department of Physics, Clarkson University, Potsdam, New York 13699-5820, USA 

ABSTRACT 

One-dimensional reaction-diffusion models A + A —»■ A + A ^ A, and A + B — > 0, 
where in the latter case like particles coagulate on encounters and move as clusters, are 
solved exactly with anisotropic hopping rates and assuming synchronous dynamics. 
Asymptotic large-time results for particle densities are derived and discussed in the 
framework of universality. 



PACS numbers: 05.40.+j, 82.20.-w 



Diffusion-limited reactions involving aggregation and annihilation processes are im- 
portant in many physical, chemical and biological phenomena [1] such as star formation, 
polymerization, recombination of charge carriers in semiconductors, soliton and anti- 
soliton annihilation, biologically competing species, etc. In this work we study by exact 
solution effects of anisotropy in ID, for single-species reactions A + A — > A or 0, and a 
two-species annihilation model A + B — > in which like particles coagulate irreversibly. 
Scaling approaches [1-2] suggest that in ID these reactions are fluctuation-dominated, 
and we cannot expect the rate equation approach to be valid. Indeed the mean-field 
rate equation approximation ignores effects of inhomogeneous fluctuations. Exact so- 
lutions and asymptotic arguments, in ID, have been used [3] to check general scaling 
and universality expectations. The ID reactions have also found some experimental 
applications [4]. These studies have assumed isotropic hopping (diffusion). 

For the reaction A + B — > 0, numerical results and phenomenological considerations 
suggest [5] that making the hopping fully directed would change the universality class in 
ID. Specifically, the large-time particle concentrations (assuming equal densities of both 
species) would scale according to c(t) ~ t -1 / 3 instead of the isotropic- hopping power 
law t -1 / 4 . A few exact and numerical results available in the literature on anisotropic 
reactions involving only one species [6] indicate that the power law is not changed. The 
model of [5] assumed that like particles interact via hard-core repulsion; this seems to 
be an essential ingredient for observing the changeover in the universality class. 

In this work we report the exact solution for two-particle annihilation with anisotropic 
hopping. However, in order to achieve exact solvability we took "sticky-particle" rather 
than hard-core interactions: the like particles coagulate on encounters and diffuse as 
groups. Our exact calculations yield the t -1 / 4 power law, found earlier for "sticky- 
particles" with different dynamics and isotropic hopping [7] . For unequal initial concen- 
trations, the large-time behavior changes; the crossover between the two regimes is de- 
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rived analytically. We also obtain exact results for A + A — > A or with anisotropy. The 
universality class of dynamics of these reactions is not affected by hopping anisotropy. 
An extended version of this work will be reported elsewhere [8]. 

In lattice models the particles hop randomly, to the extent allowed by their inter- 
actions, to their nearest neighbor sites. Two like particles can annihilate on encounters, 
A + A — > 0, or aggregate, A + A — > A. The ID kinetics of these reactions is non-mean- 
field, with the typical large time diffusional behavior of the concentration (density per 
site), c(t) ~ t -1 / 2 . For the two-species model, to be termed the AB model, unlike parti- 
cles annihilate, A + B — > 0. When like species meet, some interaction must be assumed. 
The simplest interaction is hard-core. Assuming equal A- and S-concentrations and 
random, uniform initial conditions, particle concentrations in the isotropic case scale 
according to c(t) ~ t -1 / 4 in ID. A surprising recent result [5] is the new exponent 
~ 1/3, replacing 1/4, for anisotropic hard-core particle hopping. 

In order to obtain a solvable model in ID, we consider the AB annihilation model 
with the "sticky particle" interaction. Thus, like particles coagulate irreversibly on 
encounters, e.g., nA + mA — > (n + m)A, and diffuse as clusters. When unlike clusters 
meet at a lattice site, the outcome of the reaction is nA + mB — > (n — m)A if n > m, if 
n = m, and (m — n)B if n < m. Recent numerical results and scaling considerations for 
these reactions [9] in D = 1, 2, 3 indicate that they are mean-field in D = 2, 3. However, 
in ID the power-law exponent for the density is 1/4 [7,9], with a faster power-law decay 
~ t -3 / 2 of the minority species in case of unequal densities of A and B. 

Following [10], we first consider diffusion of nonnegative charges on the ID lattice. 
Initially, at t = 0, we place positive unit charge at each site with probability p or zero 
charge with probability 1 — p. Furthermore, we consider synchronous dynamics, i.e., 
charges at all lattice sites hop simultaneously in each time step t — > t + 1, where the 
probabilities of hopping to the right, r, and to the left, i = 1 — r, are not necessarily 
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equal. This dynamics decouples the even-odd and odd-even space-time sublattices; it 
suffices to consider only those charges which are at the lattice sites j = 0, ±2, ±4, ... at 
times t = 0, 2, 4, . . ., and lattice sites j = ±1, ±3, ±5, ... at times t = 1, 3, 5, — The 
"interaction" between the charges is defined by the rule that all charge accumulated at 
site j at time t coagulates. There can be 0, 1 or 2 such charges arriving at j, depending 
on the random decisions regarding the directions of hopping from sites j ± 1. 

This model can also be viewed as diffusion-coagulation of unit-charge "particles" 
C, i.e., nC + mC — > (n + m)C. Such reactions, without the limitation of positive or 
integer charges, and with an added process of feeding-in charge at each time step, have 
been considered as models of self-organized criticality and coagulation [11-12], assuming 
isotropic hopping, r = £ = \. This coagulation reaction can be mapped [7,10] onto both 
our single-species and "sticky" AB models. However, before discussing this mapping, 
let us present the exact solution of the model of coagulating charges with anisotropic 
hopping, following the ideas of [10-11]. 

We define stochastic variables, Tj(t) = 1 or 0, with probabilities r and £, respec- 
tively. The stochastic equation of motion for the charges qj(t), equal to the number of 
C particles at site j at time t, is 

q n (t+l) = T n - 1 (t)q n - 1 (t) + [l-T n+1 (t)]q n+1 (t) . (1) 

The total number of C-particles, or the total charge, in an interval of k consecutive 
proper-parity-sublattice sites, starting at site j at time t, is given by 

fc-i 

S k,j(t) = ^qj+2i(t) = qj(t) + q J+2 {t) + ■■■ + +q l+2k - 2 {t) . (2) 

i=0 

Due to conservation of charge, the equations of motion (1) yield the following relation, 



S k ,n(t + 1) = T n _i(t)<? n _i(t)+<? n+ i(t)H h<?n+2fc-3 W + [l ~ T n+2 k-l (*)] <?n+2fc-l (t) . (3) 
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Thus, only the two random decisions at the end points are involved in the dynamics of 
charges in consecutive-site intervals. The exact solvability of coagulating- charge models 
is based on this property [11]. 

Let us introduce the function I(s, m) = S SjTn , and averages, fk,m(t) = (I (Sk,j(t),m)). 
The averaging is over the stochastic dynamics, i.e., over Tj(£), as well as over the initial 
conditions. Since the latter are uniform, fk, m (t) do not depend on j. Other choices for 
I(s, to) have been used [10-12]. In our case fk,m(t) correspond to the probability to find 
to charge units in an interval of k sites, so that /i, m (t) is the density (fraction) of sites 
with charge to. 

The variables T;(t) and Sk, n (t) are statistically independent because the latter de- 
pend only on "decision making" at earlier times. As a result, one obtains, by (3), the 
discrete diffusion-like equation, 

fk,m(t + l)=r£ [fk+l,m(t) + fk-l,m(t)] + + **) /*,»»(*) • (4) 

The m-dependence only enters via the initial conditions: fk,m(fy = p m (l — P) k ~ m {^) , 
provided < to < k, and for m > k. We also define the boundary conditions 
fo,m(t) = 1(0, to) = 5o :m in order to extend (4) to all t = 0, 1, 2, — 

In order to solve (4) we introduce the double generating function, gf-(u,w) = 

E^o Em=o/ fe . m W ttitt|m - It * s a ^ so convenient to introduce the variable a = r — £ 
directly measuring the hopping anisotropy, r = (1 + a)/2 and t = (1 — a)/2. One can 
then derive the following equations, 

g k+1 (u,w) + 2 ^^ a ^ — -g k (u,w) +g k -i(u,w) = -- — (wp + 1 - p) k , (5) 
(1 — or) u (1 — a z ) u 

with the initial and boundary conditions gk(0,w) = (wp + 1 — p) k and go(u, w) = 
1/(1 -u). 

The solution of (5) is obtained as a linear combination of the special solution 
Vt(wp+1— p) k proportional to the right-hand side, and that solution of the homogeneous 
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equation which is regular at u = 0. The coefficient Q is obtained by substitution, 

= 4(wp + l-p) 

(1 — a 2 ) it (tup + 1 — p — A_|_) (wp + 1 — p — A_) 

where A± are the roots of the characteristic equation, 

2-(l + a 2 )n±2y/(l-n) (1 - a*u) 

A± = ' {) 

The root A_, which is nonsingular as u — > 0, gives the homogeneous solution pro- 
portional to A_, where the proportionality constant is determined by the boundary 
conditions. In summary, the solution takes the form 

g k (u, w) = - fi) A* + Q(wp + l-p) k . (8) 

Densities of reactants at lattice sites derive from fk=i,m(t)- The m-dependence here 
follows by expanding (8) in powers of w. The resulting -u-dependence is complicated. 
Therefore we will keep the time-dependence in the generating-function form. Our ex- 
plicit time-dependent expressions will be derived as asymptotic results valid for large 
times. The power series in u are then controlled by the singularity at u = 1, and analyt- 
ical results can be derived by appropriate expansions. We consider the time-generating 
function for the quantities /i, m (t) which represent the probability to find charge m at 
a lattice site at time t, 

4A+(-p) m 



t=0 



1 — u (1 — a 2 ) u 



(l-a 2 )u(l-p- A+y 



We now turn to the single-species reactions introduced earlier. Our approach follows 
recent work [10] and related ideas, e.g., [13]. Consider first the reaction A + A — > A. 
In the coagulating-charge model we now regard each "charged" site as occupied by an 
A-particle, and each "uncharged" site as empty of A-particles. The dynamics of the 
coagulating charges then maps onto the dynamics of the reaction A + A — > A. The 
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quantity fi,o(t) gives the density of empty sites in both models. Therefore, the particle 
density (per lattice site), c(t), in the aggregation model, is given by c(t) = 1 — fi,o{t), 
where c(0) = p. The generating function follows from (9), 

E(u) = jr,c(t)u t = J—-Go(u) = ^^ + - * {1 ~ P) . (10) 

^—^ 1—u 1 — u (1 — or) u (1 — p — A + ) 

The function E(u) is regular at u = 0; the Taylor series is controlled by the singu- 
larity at u = 1, 



E(u) 



1 +o(i) 



(11) 



This yields the leading-order large-time behavior, 



c(t) 53 Tfr^r ' (12) 

We are not aware of other exact solutions for this model with anisotropic hopping. 
However, the leading-order large time behavior is expected to be universal in that it 
does not depend on the initial density p. Furthermore, the particle diffusion constant 
V(a) = (l — a 2 ) T>(0) decreases proportional to 1— a 2 when the anisotropy is introduced. 
Therefore, as a function of V(a)t, the result (12) does not depend on the anisotropy and 
in fact it is the same as expressions found for other A + A — > A models, with different 
detailed dynamical rules [3] . 

For the reaction A + A — > 0, the appropriate mapping is to identify odd charges 
with particles A and even charges with empty sites [10]. The generating function is 
obtained as follows, 

E(u) =jTG 2j+1 (u) = j . (13) 

j=o (l-a 2 )u (l-p-A + ) -p 2 

The large-time behavior is similar to the aggregation reaction, with the universal ex- 
pression which is less than (12) by a factor of 2. 
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The finite-time results for both models do depend on details of the dynamical rules. 
For our particular choice of synchronous dynamics, there exists an exact relation [10] 
which holds also for the anisotropic case, checked by comparing the generating functions, 

2c (t;p) =c A (t;2p) . (14) 

Here the subscripts denote the outcome of the reaction while the added argument stands 
for the initial density. 

For the AB model, we assume that initially particles are placed with density p, 
but now a fraction a of them are type A, and a fraction j3 = 1 — a are type B. The 
concentration difference is constant during the reaction; it remains (a — (3)p. At large 
times, this is also the limiting value of the density of the majority species, while the 
density of the minority species vanishes. In what follows we assume a > (3 without loss 
of generality; c(t) will refer to the density of the majority species A. 

The dynamics of the AB model can be related to that of the coagulating-charge 
model by adapting the ideas of [7] . The dynamics of the "sticky" A + B — > model can 
be viewed as coagulation. Thus, we consider the AB particles as new charges, +1 for 
A, and —1 for B. If the net charge of a coagulated cluster is positive than we regard 
it as a group of A particles (equal in their number to the charge value). If the charge 
is negative, we consider the cluster S-particle, while if the charge is 0, we regard this 
cluster as nonexistent (0). 

The probability of having an m-particle cluster in the original positive-charge-only 
model was given by /i, m (t). Each such cluster can have charge n = —m, — m+2, . . . , m— 
2, m, where we now refer to the new, ± charge definition. The key observation is that 
having a "species" label assigned to a particle at time t = is statistically independent 
of its motion and coagulation as part of clusters at later times. The density (per site) 
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of m-particle clusters with exactly n units of charge can be calculated as follows, 



m-\-n m — n 

*m,n(*) = OC 2 (3 ^ 



mi 



(15) 



/ m+n \ | / m—n \ \ 
\ 2 J' \ 2 )• 

The concentration of A-particles, i.e., the density per site of the + charge, can be written 
as c(t) = J2n°=i n J2m=n,n+2,...^rn,n(t) ■ After some algebra, we get the generating 
function, 

E(u) = - , 4A+ . r (x^-- y^-\ S(x, y) . (16) 

v 1 (l-a 2 )u(p + A + -l)\dx y dyJ K ' yj K ' 

Here we introduced the function 



oo oo 



s(x iy ) = j2Y, xn+j y ifn+2j ^- 



n=l j=0 



2x 



VI - 4xy (1 - 2x + VI - ^xy) ' 



(17) 



and the variables x = pa/ (p+ A + — 1), y = pf3 / (p + A + — 1). The evaluation of the 
double-sum is quite nontrivial; see [8] for details. 

It is useful to introduce the parameter b = a — f3 > which measures the excess of 
A at time t = 0, 

a=(l + b)/2 and (3 = (1 - b)/2 . (18) 

For the equal concentration case, 6 = 0, the large-time behavior is governed by the 
singularity at u = 1, 



E(u) 



[2(l-a 2 ) 1/4 4Vp(l-a 2 ) 3/41 ' 1 ' 



(19) 



(1-w) 3 / 4 

The leading-order behavior of the density follows from the first term, 

c(t) w ^ 1M . (20) 

2r(3/4) (l-a 2 ) 1/4 tV4 

The most significant feature of this result is that, similar to the single-species reactions, 
the anisotropy, a, dependence can be fully absorbed in the diffusion constant, in terms 
of V(a)t = (1 - a 2 ) V(0)t. The exponent 1/4 was derived in [7] for different (isotropic) 
dynamical rules. 
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An expansion for fixed b > yields 

The leading term corresponds to the constant contribution c(t) = bp + . . . which is 
expected for the majority species. In fact, expansions near u = 1 are nonuniform in the 
limits b — > + and 6 — > 0~; here we used for the first time the fact that the majority 
species is A. The approach to the constant asymptotic density is given by the third 
term, 

(1 - b 2 ) (2 - b 2 p) , 
c t - ~ — ^ • 22 

This difference is just the density of the minority species B. As before, the anisotropy 
dependence of this leading-order power-law correction is fully absorbed in the diffusion 
rate, while the exponent is consistent with the results of [7] . 

It is of interest to explore the nonuniform behavior near 6 = within the crossover 
scaling formulation. The appropriate scaling combination turns out to be proportional 
to 6/(1 — -a) 1 / 4 , as determined by inspection of various limiting expressions. It proves 
convenient to absorb certain constants into the precise definition of the scaling combi- 
nation <7, 

a = VP (1- a 2 ) 1/4 6/(1 -up* . (23) 

In the double-limit 6^0 and u — > 1 _ , taken with fixed values of a, we obtain the 
scaling relation 

E{u) w p' 1 (1 - a 2 )' 1 b~ 3 R(a) , (24) 
where the scaling function R, analytic at a = 0, can be derived exactly, 

= °*(° + ^?) 2 . ( 25 ) 

For a«l the following expansion applies, R(a) = ^a 3 + ^a 4 + (cr 5 ). The leading 
term here reproduces the first term in (19). The latter was the limiting form for u — > 1 



at 6 = 0. The second term in (19), however, is not of the form ~ 6 3 a 4 . Corrections to 
the leading scaling behavior correspond to this term in the 6 = expansion (19). 

In the opposite limit, a — > +oo, we get the expansion R(a) = a 4 + 1 — Aa~ 2 + 
O (cr -4 ). The leading term here reproduces the first term in (21); the limit a — > +oo 
corresponds to u — > 1 at fixed small positive b. Interestingly, the next two terms in 
(21) are also reproduced in their small-6 form by the next two terms here. The second 
term yields 1/ [(l — a 2 ) b 3 p\ in E(u). The third term in (21) is reproduced with the 
numerator 4 which is the small-6 limiting value. 

The scaling description provides a uniform limiting approximation in the double- 
limit 6^0 and u — > 1. Specifically, the region of nonuniform behavior near 6 = 
is exploded by the large factor ~ (1 — w) -1 / 4 . In terms of cr, the behavior is smooth 
and well defined. For instance, the result (25) applies equally well for a < which 
corresponds to A becoming the minority species. The limit of u — > 1~ at small fixed 
6 < is described by the limit a — > — oo. The appropriate expansion takes the form 
R(<j) = —1 + 4a~ 2 + O (cr -4 ), similar in structure to the a — > +oo expansion but 
without the constant-density first term. 

In summary, our exact results for the leading-order large-time particle densities of 
reaction-diffusion models in ID show expected universal power-law behaviors. Anisotropy 
of hopping has no effect on the universality class of the models studied; it can be ab- 
sorbed in the diffusion constant. Finite-time results are more sensitive [6] to the value 
of the anisotropy parameter a; they are cumbersome to derive and of less interest than 
the leading-order expressions. One interesting exception is the relation (14) which holds 
in our synchronous-dynamics models. 
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